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We compute the wall velocity in the MSSM with W, tops and stops con- 
tributing to the friction. In a wide range of parameters including those which 
Q H ' fulfil the requirements of baryogenesis we find a wall velocity of order v w ~ 10 -2 

^ ■ much below the SM value. 

CD 

1 Bubble Wall Equation of Motion in the 
MSSM 



Energy conservation leads to the equations of motion of an electroweak bubble 
wall interacting with a hot plasma of particles: 

dV T (h 1 ,h 2 ) sr^ dmj f d 3 p 

uhl+ &K X = °' (1) 

n , dVriMM , ^dm) f d 3 p 

uh2+ dh 2 + 22m^ J J^2Ei Sf ^ x) = °- (2) 

where /, = /o,i + Sfi is the distribution function for a particle species in the heat 
bath |Q,|2|. We have to sum over all particle species i. The distribution function 
is divided up into equilibrium part /o,i and out-of-equilibrium part 5 fa. The 
equilibrium part has been absorbed into the equilibrium temperature dependent 
effective potential Vr(hi, fi2). 

In the following we will restrict ourselves to late times leading to a stationarily 
moving domain wall where the friction stops the bubble wall acceleration. This 
is a reasonable assumption for the late stage of the bubble expansion where 
baryogenesis takes place. 



1 To appear in the Proceedings of SEWM2000, Marseille, June 14-17, 2000 



2 Fluid Equations 



The deviations <5/j from the equilibrium population density are originating by a 
moving wall. They are derived by Boltzmann equations in the fluid frame: 

<kfi = dtfi + x^-fi + p x ^-fi = -C[fi], (3) 
dx d Px 

with the population density fi and energy E = \/p% + m 2 (x). C[fi] repre- 
sents the scattering integral. The classical (WKB) approximation is valid for 
p 2> 1/L W ("thick wall"). For particles with E,p>gT this should be fulfilled. 
Infrared particles are supposed not to contribute to the friction. This is a crude 
approximation and there are additional contributions || which further lower 
the wall velocity. In the MSSM the wall thickness L w is of order 15/T-40/T, 
as found in and L w S> 1/T is fulfilled. Those particles which couple very 

weakly to the Higgs are denoted as "light particles" . Particles coupling strongly 
to the Higgs are heavy in the Higgs phase and therefore called "heavy", "su- 
perheavy" particles as the "left handed" stops do not appear in our calculation 
besides their remnants in the effective potential. We treat as "heavy" particles 
top quarks, (right handed) stops, and W bosons. The Higgses are left out. Fur- 
ther contributions produce an even smaller velocity. We assume now that the 
interaction between wall and particle plasma is the origin of small perturbations 
from equilibrium. We will treat perturbations in the temperature 5T, velocity 
6v and chemical potential 5fi and linearize the resulting fluid equations. Then 
the full population density fi of a particle species i in the fluid frame is given by 

U = expj^jil 

where we have generally space dependent perturbations 5% from equilibrium. In 
principle one must include perturbations for each particle species. A simplifica- 
tion is to treat all the "light" particle species as one common background fluid. 
This background fluid obtains common perturbations 5vb g in the velocity and 
5Tb g in the temperature. This leads us to 



5m + — (5Ti + 5T bg ) +p x (5vi + 5v bg ) 



(5) 



for the "heavy" particles. The spatial profiles of all these perturbations depend 
on the microscopic physics. We treat particles and antiparticles as one species 
neglecting CP violation which is a minor effect on the friction. It were, of course, 
important for the calculation of the baryon asymmetry. 

Since the perturbations are Lagrangian multipliers for particle number, en- 
ergy, and momentum, we can expand (||) to a set of three equations, called 
"fluid equations", coupled by the collision term C[8fj,,8T,gdv]. Performing the 
integrals p|,|2| leads to the general pattern 



Jj^f2C[f]=5 f iF lll +5TT Tl , J-^P-EC[f} = 5^T^+5TT T2} 

d " P Px C[f] = 5vTT v , (fi) 



(2vr) 3 T 3J 



For a stationary wall we can use dtdi — > r yi V v w S' i , and d z S{ — > 7«,^, where the 
prime denotes the derivative with respect to z = j w (x — v w t). Our equations 
are similar to those in Q but there are important additional terms. 

For each heavy particle species in the plasma we have three fluid equations. 
The final form of the fluid equations can be written in a matrix notation: 

AS' + T8 = F, (7) 

where r = Tq + I/64M. The matrices A, T, To, and M can be found in Q. The 
number C4 is the heat capacity of the plasma C4 = 78c4_ + 37c4+ including light 
quarks, leptons, and sleptons in the plasma. The perturbations are combined in 
a vector 8, the driving terms are combined in the vector F. The driving term 
containing (m 2 )' can be split up into different contributions 

, 9, , dm 2 , , dm 2 , , 

The vectors 8 and F for k particle species read (index x denotes + or — , for 
fermions and bosons, respectively, for the ith particle): 

5 = [8m 8T1 T8v\ ... 5/j k 8T k T5v k ] T , (9) 

rp _ JwVw r /•_2\/ „ f m 2\> n „ f~,2\t „ f^2\l n l T 



2T 



[ci±K)' C2±(m?)' ... ci±(m2)' C2± (m^)' 0] . (10) 



where c\±, C2± denote the fermionic(+) or bosonic(-) statistical factors, respec- 
tively, defined through c n ± = f E n ~ 2 /T n+1 (f (±)')d 3 p/(2iT) 3 . Eq. (@) has to be 
solved for 8. To a first approximation, neglecting 8', we obtain 5 = r -1 .F. Then, 
including (right-handed) stop-, top- and W particles, the equations of motion 
can be rewritten in the fluid picture as 

h 2 h 2 
h" - Vt = nilwVw—h'^ h'2 - 14 = Wio^y^" ^ 

with slightly tan/3-dependent friction constants r/i i 2=T/4Gi j 2r~ 1 Fip, with con- 
stant vectors @-Fi,2 and G\ t 2- Perhaps 8' is not negligible, so we have to solve (0) 
numerically. We compare both resulting velocities later on in Fig. 



3 Wall Velocity in the MSSM 



In order to solve eqs. (|ll| ) we derive a virial theorem, based on the necessity that 
for a stationary wall the pressure to the wall surface is balanced by the friction. 
The pressure on a free bubble wall can be obtained from l.h.s. of the equations 
of motion (11) integrated with h[ 2 , e.g. 

pi = r i k " ~ H ' idx =avt= 1 ^ v ^Y {h>i)2 - (i2) 

AVt is the difference in the effective potential values at the transition tempera- 
ture T n , which is basically the nucleation temperature. Both of the Higgs fields 
develop friction terms and we have to add the pressure contributions. In the 
MSSM the approximation |I]J| to the solution by a kink h(x) = hcntjl (1 + tanh(x/L U) )) 
is a rather good choice. With hi = sin/3 and h\ = cos/3 we are left with 

_ 20L W AV T (T n )T n 
Knt sm 4 /3(r? 2 + m cot 4 /3) 




Figure 1: Wall velocity in dependence on the parameter tan/3 for mu = —60, 0, 60 GeV; 
left: 5' = 0, right: the same plus velocities for 5' ^ resulting from the full solution 

of®. 




Figure 2: /3(x) for ml = -(QOGeV) 2 , m Q = 2TeV, m A = 400GeV, tan/? = 2, and no 
mixing. The back reaction of the friction is negligible small (left). Only artificially 
setting r]2 two orders of magnitude larger leads to sizable effects (right). 



The missing numbers for L w , h cr it, T n , and AV(T n ) can be independently deter- 
mined with methods described in @,|!|. We used the 1-loop resummed effective 
potential. The diagram is calculated for tuq = 2TeV, and uia = AOOGeV. We 
find a strong phase transition with v/T = .95 at 1-loop level for mu = — 60GeF 
at tan/3 = 2.0 and At = fi = (no stop mixing). At 2-loop level it is even much 
larger H permitting larger tan/3 for the same strength. Also with mixing in the 
stop matrix this allows to comply the experimental lower bound on the Higgs 
mass. In Figure [l] we show vw vs. tan/3 for three values of mu = —60, 0, +60GeV 
at zero mixing. Very heavy stops should decouple more and more which would 
lead to increasing wall velocities again. This behaviour is reproduced with the 
full numerical solution to (0), see Fig. |l[ Nevertheless for increasing mfj the 
used approximations become worse and corrections of order 0(m 2 /T 2 ) are im- 
portant. Our calculations are done with massless outer legs only taking into 
account changing plasma masses. 

Due to the effective potential which couples the equations of motion there 
may be back-reaction of the different friction contributions to 8(3/ dz leading to 
a change in A/? = max(d(3/dz). A larger A/3 were highly welcome to obtain a 
larger baryon asymmetry. This becomes even more important since we realized 
that in the MSSM transitional CP violation does not occur. Therefore we 
must exploit the explicit phases which may nevertheless be strongly restricted 
by experimental bounds. The determination of A/3 can be done numerically by 
solving eqs. (11) with extensions of known methods ||. But only for artificially 
large friction we obtain sizable effects (see Fig. ^). 
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